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Abstract 

The composite operator TT, obtained from the components of the energy-momentum tensor, 
enjoys a quite general characterization in two-dimensional quantum field theory also away 
from criticality. We use the form factor bootstrap supplemented by asymptotic conditions to 
determine its matrix elements in the sinh-Gordon model. The results extend to the breather 
sector of the sine-Gordon model and to the minimal models M.2/(2N+3) perturbed by the 
operator 4>\^. 



1 Introduction 



The ability of a quantum field theory to describe a system with infinitely many degrees of freedom 
is reflected by an infinite-dimensional operator space. In two dimensions, the detailed structure 
of the operator space at a generic fixed point of the renormalization group was revealed by the 
solution of conformal field theories pQ. It is divided into different operator families, each one 
consisting of a primary and infinitely many descendants. Within an operator family, the scaling 
dimensions differ from that of the primary by integer numbers that label different 'levels'. 

Perturbative arguments lead to the conclusion that this same structure is maintained when 
conformal invariance is broken by a perturbation producing a mass scale j2]. If the massive 
theory is integrable, the operator space can be studied non-perturbatively within the form factor 
bootstrap approach [31 It was shown for several models that the global counting of solutions 
of the form factor equations matches that expected from conformal field theory [3] . 

As for the correspondence between solutions of the form factor equations and operators, 
asymptotic conditions at high energies play a crucial role. While primary operators are naturally 
associated to the solutions with the mildest asymptotic behavior, we argued in [Hj that specific 
asymptotic conditions selecting the solutions according to the level can be identified. These were 
used in [Jj to show the isomorphism between the critical and off-critical operator spaces in the 
Lee- Yang model, level by level up to level 7. 

Asymptotic conditions, however, cannot determine completely a descendant operator in a 
massive theory. Indeed, they leave unconstrained terms which are subleading at high energies 
and depend on the way the operator is defined away from criticality. The operator TT, obtained 
from the components of the energy-momentum tensor, appears as the natural starting point in 
relation to the problem of the off-critical continuation of descendant operators. Indeed, being 
the lowest non-trivial scalar descendant of the identity, this operator allows for a quite general 
characterization in two-dimensional quantum field theory. A. Zamolodchikov showed how to 
define it away from criticality subtracting the divergences which arise in the operator product 
expansion of T and T |S]. We showed in [Hj for the massive Lee- Yang model that, with this 
information, the form factor programme outlined above allows to uniquely determine TT up 
to an additive derivative ambiguity which is intrinsic to this operator. Our results have been 
successfully compared with conformal perturbation theory in 

In this paper we address the problem of determining TT in the sinh-Gordon model. The 
essential difference with respect to the Lee- Yang case is that, while the latter is a minimal model 
with the smallest operator content (two operator families), the sinh-Gordon model possesses a 
continuous spectrum of primary operators, a circumstance that seriously complicates the iden- 
tification of specific solutions of the form factor equations. The massive Lee- Yang model is the 
first in the infinite series of the (^3-perturbed minimal models 2/ (2^+3)7 each one contain- 
ing N + 1 operator families. Due to a well known reduction mechanism jl()l I12[ these 
massive minimal models have to be recovered from sinh-Gordon under analytic continuation to 
specific imaginary values of the coupling. The form factor solution for the operator TT of the 
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sinh-Gordon model that we construct satisfies this requirement. 

The paper is organized as follows. In the next section we recall a number of facts about 
bosonic theories in two dimensions. The form factor solutions for the primary operators in 
the sinh-Gordon model are reviewed in section 3, while the solution for TT is constructed in 
section 4. Few final remarks are contained in section 5. Six appendices conclude the paper. 



2 Bosonic field with a charge at infinity 

If B\^ u is a tensor antisymmetric in the first two indices, the energy-momentum tensor T^ u of 
a quantum field theory can be modified into 

f llu = T lJiU + d x Bx fl u (2-1) 

preserving the conservation d^T^ v = and the total energy- momentum P v = f da^T^ v . For a 
neutral two-dimensional boson with action 



A= I d 2 x 



(2.2) 



the choice 1 



leads to 



B XiXV = -^Le X)1 e vp d^ (2.3) 

V 27T 



d X B Xflu = -^L (8^d u - Vlu ,n)<p (2.4) 
V 2ir 



and to the variation 



e = e + ^./fi> (2 ' 5) 



in the trace of the energy-momentum tensor 



TT 



©=2 T M- ( 2 - 6 ) 

The canonical definition = — 2V and the equation of motion Dip = dV /dip give the classical 
result 

iQ d 



e --'l- 1 + vS<*W- (2 ' 7) 



The parameter Q is dimensionless and goes under the name of "background charge" or 
"charge at infinity". It does not change the particle dynamics 2 but, inducing a modification of 
the energy-momentum tensor, essentially affects the scaling properties of the theory. 

Free massless case. In the free massless case corresponding to the action 



\ J d\ (d,pf (2.8) 



1- We denote by r/^ the flat metric tensor and by £a m the unit antisymmetric tensor in two dimensions. 
2 In particular, Q does not enter the perturbative calculations based on the Lagrangian. 
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a nonvanishing Q leaves the energy- momentum tensor traceless (O c / = G = V = 0) and the 
theory conformally invariant. The central charge is 

C = 1 - 6Q 2 (2.9) 

and the scaling dimension of the primary operators 

V a {x) = e^^W (2.10) 

is 

X a = 2a(a - Q) . (2.11) 
A derivation of these results within the formalism of this paper is given in appendix A. 

Minimal models. For real values of the background charge the bosonic model can be used 
to reproduce the minimal models of conformal field theory with central charge smaller than 1 
|14j . Indeed, the requirement that the 2/c-point conformal correlator of an operator V a ~ Vq_ q 
is nonvanishing for any positive k selects the values 

ot = a m>n = 2 [C 1 _ m)a + + (1 - n)a_] , m,n = 1,2,... (2.12) 

with 



2 

Equations (|2.9|) and (|2.11|) then reproduce the central charge 



a± = ^ • (2.13) 



C p/p , = 1-6^—^ (2.14) 

and the scaling dimensions 

(p'm - pn) 2 - (p - p') 2 

-A-m.n — ^ ; l z - iC) J 

2pp' 

of the primary operators </> m ,n in the minimal models M p / p i ^ through the identification 

a± = ±(l^\ (2.16) 



Since a+a_ = — 1, one obtains the correspondence 



-Vn,n ~ v a m n 



Va m , n = exp 



I 



(m — 1) (n — l)a_ 



<P> ■ (2-17) 



It is well known that, although genuine minimal models (i.e. those possessing a finite number 
of conformal families which form an operator space closed under operator product expansion) 
correspond to rational values of p/p' , the above formulas in fact apply to the degenerate operators 
of conformal field theory for continuous values of central charge smaller than 1. 
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Liouville theory. A deformation of Q2.8|) which does not introduce any dimensional parameter 
is obtained adding an operator which is marginal in the renormalization group sense, namely 
has X a = 2. This requirement selects V a _ ~ Vq-o. 



V a+ . To be definite we take 



Al = J d 2 x I(^) 2 + / x e ^^ 



where we defined 



(2.18) 



(2.19) 



and fi is a coupling constant. For real values of b this is the action of Liouville field theory, 
which is conformal and has been extensively studied in the literature (see e.g. |15j for a list of 
references). Notice that the condition C ; = gives the value Q c i = —i/b, which coincides with 
the exact result 

+ Q_ = -i [HtI (2.20) 



aj 



in the classical limit b — ► 0. The central charge and scaling dimensions of exponential operators 



in Liouville field theory are given by 1)2. 9JI and 1)2.1111 with Q given by 1)2.2 



Sinh-Gordon model. The sinh-Gordon model is defined by the action 
A s hG - 



d 2 x 



(2.21) 



which can be regarded as a perturbed conformal field theory in two different ways. 

The first one consists in seeing it as a deformation of the Gaussian fixed point, i.e. the 
conformal theory with C = 1. This amounts to setting Q = keeping b as a free parameter. In 
such a case, both the exponentials appearing in the action have scaling dimension —2b 2 and are 
never marginal for real values of b. They play a symmetric role and the theory is invariant 3 under 
the transformation ip — ► —ip. The trace of the energy-momentum tensor, being proportional to 
the operator which breaks conformal invariance, is = ~ fj, cosh y/87rbip. 

The second point of view consists in looking at ()2.21|) as the perturbation of the Liouville 
conformal theory 1)2.18)1 by the operator e -v^ b f scaling dimension Xn, = —2{2b 2 + 1). 
Since this is now the operator which breaks conformal invariance, we have ~ 
result which agrees with the classical expectation (f2.7|) once one uses Q c \ for Q. 



[I e 



Sine-Gordon model and its reductions. The sine-Gordon action 



A 



sG 



d 2 x 



— {d u tp) 2 — 2fx cos V8tt (3 (p 



can be obtained from (|2.21|) taking /i = // and 

P = -ib . 



(2.22) 



(2.23) 



For real values of (3 the only direct interpretation of this action as a perturbed conformal 
field theory is as a deformation of the C = 1 conformal theory through the operators Vp and V-p 
3 In this case the couplings fi and \i have the same dimension and can be made equal shifting the field. 
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with scaling dimension 2(3 2 . The perturbation is relevant and the theory is massive for f3 2 < 1. 
In this range the sine-Gordon model is known to be integrable and its factorized S-matrix is 
known exactly [TJj|. The particle spectrum consists of the soliton A and antisoliton A and, in 
the attractive range < f3 2 < 1/2, of their neutral bound states, the breathers B n with 

1 < n< Int (^j (2.24) 

and masses 

m n = 2M sin -y ; (2.25) 
here Int(x) denotes the integer part of x, M is the mass of the soliton and 

Taking Q / and looking at the sine-Gordon model as a perturbation of the conformal 
theory (|2.18|) is problematic because the action Q2.18JI becomes complex when b is imaginary. 
Formally, however, this point of view leads, through the identity 

/? = «_ = -/!, (2.27) 
V V 

to the conformal field theories with central charge Q2.14JI perturbed by the operator V-p ~ $1,3, 
namely to the action 

A Mp/p , +\ J cPx^ix) . (2.28) 

This (/^^-perturbation of the C < 1 conformal field theories is known to be integrable for any 
value of A [2] and massive for a suitable choice of the sign of A |17j . This choice is implied in 



(J2J2SJ). 

The relation between the sine-Gordon model and the action (|2.28l) suggested by these formal 
reasonings can be confirmed and put on firmer grounds within a framework known as quantum 
group reduction |10 [lll[IT2| llHj. This relies on the fact that the sine-Gordon S'-matrix commutes 
with the generators of the affine quantum group SL(2) q with q = exp(z7r//3 2 ), and that for 
rational values of j3 2 = p/p' a restriction can be operated in the space of particle states and 
operators of the model which is consistent with this algebraic structure and preserves locality. 
The quantum field theories obtained through this reduction mechanism indeed coincide with the 



perturbed minimal models ([2.28ft . 

While soliton and antisoliton transform as a doublet under the action of the quantum group, 
the breathers are scalars. This is why the reduction takes its simplest form when the space of 
states can be restricted to the breather sector. Let us recall that the amplitude for the scattering 
between the breathers B m and B n in the sine-Gordon model is 4 IK 



min(m,n)— 1 

s mn (0) = t {m+n) ^e)t lm _ n ^(e) J] (2 - 29) 

7 = 1 



4 The rapidity variables 6% parameterize energy and momentum of a particle as {pl,p\) = (m cosh 9i, msinh#i) 
m being the mass. The scattering amplitudes depend on the rapidity difference between the colliding particles. 
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where 

, tanh \(9 + ina) 

t a (9) = ^ '- . (2.30) 

tanh g {0 — ma) 

While the double poles are associated to multiscattering processes ^JE3j the simple poles 
located at 9 = i(m + ra)£/2 and 9 = i(ir — \m — n\£/2) correspond to the bound states B n+m and 
B\m-n\i respectively, propagating in the scattering channel B m B n . To be more precise, the first 
class of simple poles can be associated to the bound states B m+n only for values of £ such that 
m + n < Int(7r/£). Indeed, (|2.24j) shows that outside this range the particle B m+n is not in the 
spectrum of the model in spite of the fact that the pole in the amplitude (|2.29|) may still lie in 
the physical strip Im# S (0,7r). In this case, however, this pole can be explained in terms of a 
multiscattering process involving solitons as intermediate states. This is why, for generic values 
of £, the breather sector of the sine-Gordon model is not a self-contained bootstrap system. 
The situation becomes different when £ takes the special values 

2ir 

^ = ^—^, N = l,2,..., . (2.31) 

In this case (|2.24jl . (|2,25j) and (|2,29|) show that there are N breathers and that the formal 
identities 

m n = m 2 N+i-n (2.32) 

(9) (2.33) 

hold, so that the pole discussed above can be associated to B m+n for m + n < N and to 
£>2jv+i-m-n for Til + n > N, without any need to resort to the solitons. Hence, for the values 
(|2,31|) of the coupling, the breather sector of the sine-Gordon model provides alone a self- 
consistent factorized scattering theory and, consequently, defines an infinite series of massive 
integrable models labeled by the positive integer N. Equations (|2.26j) and (|2.27|) then identify 
these massive models with the minimal models A^2/(2Af+3) with central charge 

C N = l-^ 2 l +1>2 . P-34) 
2iV + 3 ' v ; 

perturbed by the operator <j)i 3 with scaling dimension 

l|? = -2^. (2.35) 
1,3 2iV + 3 y ' 

This conclusion was first reached in (11] . The thermodynamic Bethe ansatz |20( I3flj confirms 
that the scattering theory (|2.29f) gives the central charges (|2.34|) for £ = £at, a result that can be 
regarded as a non-perturbative confirmation of the fact that the charge at infinity does not affect 
the dynamics of the particles. The minimal models A^2/(2Af+3) possess the iV non-trivial primary 
fields <j)\ t ki k = 2, . . . , N + 1, plus the identity <j>\ % \. The negative values of the conformal data 
(|2.34l) and (|2.35|) show that these models do not satisfy reflection positivity. The case N = 1 
corresponds to the Lee- Yang model |221 1231 the simplest interacting quantum field theory. 
Its 5-matrix Sn(9)\^ =27r / 3 = t 2 /3(9) was identified in [23] . 
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3 Primary operators in the sinh- Gordon model 



Most of the results discussed for the sine-Gordon model apply to the sinh-Gordon model (|2.21|) 
through the correspondence (|2.23[) . In particular also the latter model is a massive integrable 
quantum field theory. The particle B interpolated by the scalar field corresponds to the sine- 
Gordon lightest breather B\. The scattering amplitude 



S{9) = S n {9)=t L {9) 



(3.1) 



does not possess poles in the physical strip when b is real and completely specifies the S-matrix 
of the sinh-Gordon model. This amplitude was proposed and checked against perturbation 
theory in b in [2f)J El HE] • It should be clear from the discussion of the previous section that the 
S'-matrix is the same for the two ultraviolet limits (Gaussian fixed point and Liouville theory) 
compatible with the action (|2.21j) . 

The S'-matrix determines the basic equations satisfied by the matrix elements of a local 
operator Q(x) on the asymptotic multiparticle states The form factors 5 



71, . . . ,u n 



0„) = <O|$(O)|£(0i)...B(0 n )> 



obey the equations 



F*(0 1 + a,... ) n + a 



F n . . . ,9i, 0i+i, . . . , 9 n ) — S(9i — 0«+i) F n 
F*(0i + 2i7r,0 2 ,...,0 n ) = F*(0 2 ,...,l 



H+li "ii ■ ■ ■ i "nj 



>m u l, 



Rest 



i-Y[s(e-dj) 



F*(0i, . . . 



(3.2) 



(3.3) 
(3.4) 
(3.5) 

(3.6) 



where the euclidean spin s$ is the only operator-dependent information. 

The solutions of the equations (|3.3j) - (|3.6j) can be parameterized as [29^ l2T] 



F*(0i, . . . 



u?(e 1} 



n 



i<j COsh "V* 



(3.7) 



Here the factors in the denominator introduce the annihilation poles prescribed by ()3.6j) . and 

(t7T - 0)t 



with 



^•(0) = AA(0 exp 

= - 

Af(£) = ^(m) = exp 



(it . sinh | 



2tt 



4 sinh — sinh 
2?r 



1 + 



TT 



dt . . 

7^ sinh 2 * 



sinh | 1 



(3.8) 

(3.9) 
(3.10) 



3 We denote by |0) the vacuum state. 
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The function J-{9) is the solution of the equations 

JF(0) = S(0)F(-0) 



(3.11) 



T(9 + 2m) = F{-9) (3.12) 

with asymptotic behavior 



lim F{9) = 1 ; (3.13) 



it also satisfies the functional relation 



sin Ti n 

^ + ^ = ___ . (3.14) 

The expression (|3.8j) is convergent in the range — 7r < £ < which is relevant for the sinh-Gordon 
model. 

All the information about the operator is contained in the functions £/*. They must be 
entire functions of the rapidities, symmetric and (up to a factor (— l) n_1 ) 2-7rz-periodic in all 9j. 
We write them in the form 

/ \(«-l)/2 

U*(9 1 ,..,9 n )=H n l-^j Qt(0i,.-,0 n ) (3.15) 

using the symmetric polynomials generated by 

J 

i=l k=0 
with Xi = e 9i , and choosing the constants 



f[(x + Xl ) = Y, x n - k a^ fa, ... , x n ) (3.16) 



(3.17) 



The equations (|3.3|) - ()3.6|) imply 



Qt{9 1 + a,..,e n + a) = e\ s * + 2 ) a Q* (9i,..,9 n ) (3.18) 

Qn (^1) -j 9i,9i + i, .., n ) = .., 9i + i,9i, .., n ) (3.19) 

Q* (01 + 2vri, .., n ) = Qt (9 U .., n ) (3.20) 

Q* +2 (0 + i7r,Mi,",0n) = (-l) n xD n (x,a;i,..,x n )Q*(ei,..,e n ) , (3.21) 

where x = e e and 

D B (x,x ll ..,x B ) = X; E (-l) fc+1 H^ 2(n - fe)+m 4 n) 4-L> (3-22) 

fc=l m=l,odd 

with 

r , sin(m£) 

H = -^f- (3.23) 



S 



The equations (j3.18j) - (|3.21|) admit infinitely many solutions which account for the infinitely 
many operators with spin s$. The scalar (s<j> = 0) solutions with the mildest asymptotic behavior 
are expected to correspond to the primary operators of the theory, namely the exponential 
operators (|2.1U|) . Introducing the notation 



the asymptotic factorization condition 



lim F*° (9i + X,...,9k + X, 9 k+1 , ...,6 n ) = if 9 k )F*° k (9 k+1 , ...,9 n ) (3.25) 

characterizes the scalar primary operators <3?o with non-vanishing matrix elements on any number 
of particles and in this case selects the solutions |3L)j 

Q^(9 1 ,...,9 n ) = [a]detM^(a), (3.26) 

where M^ n \a) is the (n — 1) x (n — 1) matrix with entries 

< ) (a) = [a + i-j>W (3.27) 

and a is a complex parameter. The trigonometric identity [a] 2 — [a — 1] [a + 1] = 1 is useful 
to check that the functions (|3.26j) solve the recursive equation (|3.21|) . It was found in (HO] (see 
also 32 ) that the solution (|3.26f) corresponds to the exponential operator 6 V ia b = VL Q/ g. The 
solutions ()3.26j) satisfy the property 

Q^ = (-l) n Ql a \ (3.28) 

which is expected since a — > —a amounts to (p — > —ip, and the particle B is odd under this 
transformation. In presence of the charge at infinity Q = ibtr/^ we expect the identification 
V a = VQ- a , namely 

Q { n ] = Q ( n^ K) ■ (3-29) 

The identity 

[a] =-[a + 7r/£], (3.30) 

together with ()3.28|) . ensures that (|3.29|) holds. 

The solutions for the exponential operators allow to determine the form factors of the com- 
ponents of the energy-momentum tensor. Let us start with the trace (|2.fi|) . As seen in the 
previous section, when the ultraviolet limit is taken to be the Liouville theory (|2,18|) we have 7 
© ~ Vib, so that 

6 Since the identity [a + 2tt/£] — [a] implies the periodicity Qi^ 2 "^^ = Q^n \ we consider values of a in the 
range (0, 2-7r/£). 

7 From now on we omit the tilde on the components of the energy-momentum tensor. 
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where the normalization is fixed by the condition 



F%(9 + m,9) = |m 2 (3.32) 

which corresponds to the normalization 

(B(8)\B(e')) = 2vr 6(9 - 6') (3.33) 

of the asymptotic states. If instead the sinh-Gordon model is seen as a perturbation of the 
Gaussian fixed point, the symmetric combination Qn + Qn ^ must be considered. Then (|3.28f) 
shows that the result 1)3.31(1 still holds for n even, while Q® vanishes for n odd. This analysis 
was first performed in [SHI, where the ultraviolet central charge was also evaluated through the 
C-theorem sum rule ^3 [HI] 

C = H J d 2 x \x\ 2 (e(x)@(0)) conn . (3.34) 

It was checked that a truncated spectral expansion of the two-point trace correlator in terms of 
the form factors (|3.31|) reproduces with good approximation the Liouville central charge (|2.9|) if 
the sum is performed over all n, and the Gaussian value 1 if the sum is restricted to the even 
contributions [3*3] . 

When inserted in the asymptotic factorization equation (|3.25(l with = the solution (|3.31|) 
prescribes the result 

<e>-tf~J^, ,3.35) 

which coincides with that known from the thermodynamic Bethe ansatz (see |35j ) . 

The form factors of the other components of the energy-momentum tensor are easily obtained 
exploiting the conservation equations 8 

BT = dG 

Bf = BQ (3.36) 



which lead to 



a (n) a (n) q 

F n . . . , 9 n ) = 7-7 — F n ( 



(7 



(n) 
n-l 

An) 



FZid!, ...,9 n ) = („" \nj F n(.8li ■■-,&n) (3.37) 

for n > 0; (T) = (T) = as for any operator with nonzero spin. 

It follows from the discussion of the previous section that, through the analytic continuation 
(|2.23|) . the above results for the exponential operators in the sinh-Gordon model also hold for 
the matrix elements of these operators on the breather B\ of the sine-Gordon model. Moreover, 
8 We use the notation d = d z and 5 = dz with reference to the complex coordinates z = Xi+ix2 and z — x\ ~ixi- 
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when the coupling £ takes the discrete values H'2.31j) corresponding to the reduction to the </>i,3- 
perturbed minimal models A^2/(2Ar+3)> these results give, through the correspondence (|2.17j) . 
the form factors for the independent primary operators (f>n, I = 1, . . . , N + 1, of these massive 
minimal models. The reduction from the continuous spectrum of exponential operators of the 
sinh-Gordon model to the finite discrete spectrum of primary fields in the massive minimal 
models follows from the fact that in the latter case the form factors have to satisfy constraints 
on the bound state poles in addition to (|3.3j) - (|3.6|) . For example, the fusion B\B\ — > B2 requires 

Res e > = e +ii F^. 2 (9',e,e 1 , . . . ,0„) = ir?i(O|*(O)|B 2 (0") B i(0i) • • • B M) , (3-38) 
where 6" = 6 + i£/2 and the three-particle coupling is obtained from 

Res e= i 4 Sn(0) = i {rjj 2 . (3.39) 

It turns out ^2 EE] that for £ = £at the complete set of bound state equations implied by the 
breather sector seen as a self-contained bootstrap system selects among the solutions 9 (EOfl) 
only those with a = 1, . . . , N, besides the identity (see appendix[E]). These solutions correspond 
to the operators 

fc-fcft* = &,2k+l , k=l,...,N (3.40) 
{13 N = -y/2/(2N + 3)), which, in view of the reflection relation 

h,l = h,2N+z-u l = l,...,2N + 2 (3.41) 

are all the primaries of this series of minimal models. The results (|3.31|) . I|3.35() and (|3.37|) for 
the matrix elements of the energy-momentum tensor apply to the minimal massive models for 

4 The operator TT 

We have seen in the previous section how the primary operators correspond to the solutions of 
the form factor equations (|3.3|) - (|3.6|) with the mildest asymptotic behavior at high energies. The 
remaining solutions of these equations should span the space of descendant operators. 

At criticality descendant operators are obtained acting on a primary with products of Vira- 
soro generators L_j and L-j- The sum of the positive integers i (j) defines the right (left) level 
I (I) of the descendant. We denote by a descendant of level (I, I) of a primary $0,0 = ^o- 
Due to the isomorphism between critical and off-critical operator spaces, the notion of level 
holds also in the massive theory. A relation between the level and the asymptotic behavior of 
form factors has been introduced in El • in particular, for operators $2^ with non-zero matrix 
elements on any number of particles this relation reads 

Fn 1 ' 1 (fli + A, ... ,fl fc + A, 6 k+1 , ...,6 n )~e l \ A — > +00 (4.1) 

9 We stress that for £ = £at the knowledge of the Q„ for all n completely determines the operator since all 
matrix elements involving particles Bj with j > 1 can be obtained through fusion equations like |3.88p . 
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for n > 1 and 1 < k < n — 1. It was also argued that the scaling operators of this kind satisfy 
the asymptotic factorization property 

lim e~ lx F£ (0i + A, ... A + A, fc+1> . . . , n ) = j^F^ (0 X , . . . , fc ) 0?k+i. ■ • • , 5 n ) , 

(4.2) 

where and Ci are operators that in the conformal limit converge to the product of right and 
left Virasoro generators, respectively, acting on the primary. This equation reduces to (|3.25l) in 
the case of primary operators. 

This asymptotic information can be used to classify the form factor solutions according to 
the level and to determine the leading part in the high-energy (conformal) limit. The subleading 
contributions depend instead on the way the operators are defined off-criticality. A complete 
analysis of the operator space in the massive Lee- Yang model up to level 7 is given in [7] . 

The composite operator TT, obtained from the non-scalar components of the energy-momentum 
tensor, is the simplest non-derivative scalar descendant of the identity. A. Zamolodchikov showed 
that this operator can be defined away from criticality as [S] 

TT(x) = lim[T(x 4- e)T(x) — @(x + e)0(x) + derivative terms] , (4-3) 

where 'derivative terms' means terms containing powers of e times local operators which are total 
derivatives. One consequence of this equation is that, if \n) and \m) denote n- and m-particle 
states with the same energy (E n = E m ) and momentum (P n = P m ), the equation 

(m\TT(0)\n) = (m\T(x)T(0)\n) - (m|0(ar)0(O)|n) (4.4) 

holds, with the r.h.s. that does not depend on x. Since the generic matrix element can be 
reduced to the form factors ()3.2|) iterating the crossing relation 

(B(9' m ) . . . B(e[)\m\B(9i) ■ ■ ■ B(9 n )) = (B(9' m ) . . . B(9' 2 )\m\B(9[ + m)B{9 1 ) . . . B{9 n )) + 

n i— 1 

2vr^ S(9[ - 6$ [J S(9 k - 9[) (B(9' m ) . . . fl(0£)|*(O)|B(0i) . . . B(9^)B(9 l+1 ) . . . B(9 n )) , 

i=l k=l 

(4.5) 

the identities (|4.4|) contribute to the identification of the form factor solution for the operator TT, 
in particular of the subleading parts which are left unconstrained by the asymptotic factorization 
property Q4.2|) . Since TT is a level (2,2) descendant of the identity, the factorization takes in 
this case the form 

lim e- 2X F? f (0i + A, . . . , 9 k + A, 9 k+1 , . . . , 9 n ) = F fc T (0i, . . . , 9 k ) F*_ k (9 k+1 , . . . , 9 n ) . (4.6) 

We now argue that ()4.4I) implies the property 

F^ = sum of terms containing (cr^yia^ 1 }^ , i I + j > 2 (4-7) 

for n larger than 2 but other than 4. Observe first that the components of the energy-momentum 
tensor, being local operators of the theory, must have form factors whose only singularities in 
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rapidity space are the annihilation poles prescribed by (|3.6[) . and possible bound state poles 37 . 
Then, it follows from 1)3. 3 7 j) that 

Fjoc(a^) 2 , Fjoc^) 2 , ^oce^o^, (4.8) 

for n > 2. Use now the resolution of the identity 

k=0 J 

to expand the r.h.s. of ()4.4|) over matrix elements of T, T and @. If the total energy- momentum 
of the intermediate state \k) differs from that of \m) (which, we recall, equals that of |re)) the 
two sums in the r.h.s. separately depend on x, and must cancel each other in order to ensure the 
^-independence of the result. Then we are left with the contributions of matrix elements over 
states with the same energy and momentum, which all are x-independent. Consider the case m 
and n both larger than zero, to ^ n. It follows from (|4.8j) that each of these matrix elements 
generically vanishes at least as r/ 2 , if 77 is an infinitesimal splitting between the energies of the 
two states in the matrix element. A milder behavior as 77 — > is obtained when \k) is identical 
to I to) or to \n). Indeed it can be shown using (|4.5|) and (J3.6|) that the matrix elements of the 
energy-momentum tensor over identical states are finite and non-zero. In the r.h.s. of (|4.4|) 
these non-zero matrix elements multiply a matrix element vanishing at least as rj 2 . We argue 
in a moment that the form of the 'derivative terms' in (|4.3|) is such that they contribute terms 
vanishing at least as rj 2 to the r.h.s. of (|4,4j) . Putting all together, we conclude that the l.h.s. 
vanishes at least as ry 2 , and this implies the form (|4.7|) . A more detailed derivation including the 
explanation of the limitations on n can be found in appendix B together with the form factor 
expansion of (|4.4I) . 

The property (|4.7f) can also be understood in the following way. Since mcrj and ma^ l } 1 /a^ 
are the eigenvalues of P = id and P = —id on an n-particle asymptotic state, 1)4.8)1 follows 
from the fact that in two dimensions the energy-momentum tensor can formally be written as 
(2/7r)(r? Atv D - d^d u )A(x), or 

T = d 2 A , f = B 2 A, Q = ddA, (4.10) 

in terms of an operator A(x) which is not a local operator of the theory 10 . Using the notation 
A ■ B = A{x + e)B(x) one has 

T ■ f - 6 • 9 = \d 2 {A ■ f) + \d 2 (A ■ T) - dd(A ■ 6) . (4.11) 

The property (|4,7|) then follows from the fact that also the 'derivative terms' in (|4.3|) must be 

derivative operators of at least second order: those associated to negative powers of e because 

they must cancel the divergences arising in (|4.11|) when e — > 0; those associated to e° because 

they must have I = I <2 and we know that L_i = d, L—\ = d. 

10 Indeed, F^(9i, 62) contains a pole at 61—62= ivr, in contrast with 113.611 which prescribes a vanishing residue. 
Essentially, this is why l|4.8^ holds only for n > 2. 
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Having collected this information, we can move forward in the determination of the form 
factors of TT. The first requirement to be satisfied is that they solve the form factor equations 
(|3,3|) - (|3,6|) and have the asymptotic behavior (|4.1|) with I = 2. Since the equations ()3.3|) - (|3,6|) 
are linear in the operator, the solution for TT can be written as a linear superposition of the 
scalar form factor solutions behaving as in (|4,1[) with I = 0, 1, 2. We will refer to these solutions 
as level 0, level 1 and level 2 solutions, respectively. We now argue that this linear combination 
can be restricted to 

Fl f = am- 2 Ff B ' e + if + cFf e + dm 2 F* + em 4 F n 7 , (4.12) 

with a, c, d, e dimensionless constants and F^ to be defined below. Indeed, the level solutions 
corresponding to scaling operators are spanned by the primaries and, among these, only the 
form factors 

Fn = Ko (4-13) 

of the identity and those of G satisfy the requirement Q4.7JI . As for the level 1 scaling operators, 
they are all of the form L_iZ_i<l>o = dd&o. Equations Q4.4|) . (|4.6|) and Q4.7JI put no constraint 
on the contribution to (|4.12j) of such level 1 derivative operators. However, we expect that the 
form factors of the operator TT in the sinh-Gordon model enjoy the following properties 11 : 

i) to give the form factors of TT on the lightest breather of the sine-Gordon model under 
analytic continuation to positive values of £; 

ii) to give the form factors of TT on the lightest particle of the </>i ^-perturbed minimal 
models .^2/(2^+3) when we set £ to the values £n defined by (|2.31l) : 

iii) to be continuous functions of ^. 

We know that, while the sinh-Gordon and sine-Gordon models possess a continuous spectrum of 
primary operators, the ^i^-perturbed minimal models A^2/(2iV+3) possess a discrete spectrum 
of N + 1 primaries. The identity and the trace of the energy-momentum tensor are the only 
primary operators which are present in all these models. Then the natural way to comply with 
the requirements i)— iii) is that I and are the only primaries that contribute to (|4.12fl . The 
extension of the argument to any of the models (|2.28|) with p/p' rational implies that this is 
actually the only possibility. While we had reached this conclusion about the primaries by 
another path, the present reasoning also requires that dd® and d 2 d 2 @ are the only derivative 
operators which can appear in (|4,12|) , 

The superposition (|4.12|) without the term F^ is sufficient to provide the most general 
parameterization with the required asymptotic behavior up to n = 2. This is why F^ is a scalar 
three-particle kernel solution of the form factor equations, namely has the property 

if = 0, n = 0,1,2; (4.14) 

the first two non- vanishing elements of this solution are F^ and Fjf, with and which 
factorize Tli<j cos h ' 2 3 i n such a way to satisfy (|3,6[) with on the r.h.s. The F£ are made of 
11 Analogous properties do hold for the components T, T and Q of the energy-momentum tensor. 
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terms which under the limit (|4.1[) behave as e with I = 0, 1, 2. They define the local operator 
K,{x) which in (|4.12j) accounts for the linear independence of TT within the operator space of 
the theory. 

Since our equations do not constrain the level 1 derivative contributions to (|4.12|) . the coef- 
ficient c will remain undetermined. This conclusion agrees with conformal pertubation theory, 
which states that, due to the resonance phenomenon [2] (see also |38j ) . the operator TT can 
only be defined up to a term proportional to ddO (x) [HI- Our remaining task is that of showing 
that the r.h.s. of (|4.12|) can be uniquely determined up to this ambiguity. 

For n = 2, the second derivative term is the only one contributing to the limit in (|4.6|) . and 
this fixes 

The coefficients d and e can be determined by (|4.4[) with m = n = 1. Indeed, when we use 
(J4.9J) to expand the operator product in the r.h.s., the x-independence of the result implies that 
only the k = 1 intermediate state gives a non-vanishing contribution. Using (|4.5j) to go to form 
factors we obtain the identities 12 

F 2 Tf (*7T, 0) = -2(G)F 2 °(i7r, 0) = -7rm 2 (6) (4.16) 



(TT) = -(G) 2 . (4.17) 
On the other hand, we have from (|4,12|) that 

(TT) =dm 2 (Q) +em A (4.18) 



F 2 TT (i7r,0) = dm 2 F 2 (ivr,O) , (4.19) 

so that we obtain 

d = -- ^ (6) (4.20) 

-5- (4.2!) 

The search for the kernel contribution to (|4.12j) starts from the most general solution 

of the form factor equations satisfying (|4.14j) . 1)4. 7 Jl and (|4.1|) with I = 2. We checked that 

equations (|4.6|) and (|4.4I) uniquely fix the level 2 and level parts, respectively, within such a 

solution. Concerning the level 1 part, it cannot be determined by these conditions, because they 

do not exclude the contribution of those linear combinations of <9<9-derivatives of the primaries 

which vanish on one- and two-particle states. This indetermination is eliminated if we impose 

the conditions i)-iii) above. To see this consider the c^i^-perturbed minimal models M2/(2N+3)- 

12 The identity l|4.17|l was originally observed in |.S8) and follows also from l|4.4^ with m — n — and \x\ — > oo. 
See also |39] for results on the vacuum expectation values of descendant operators in integrable models. 
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The identification B2N+i—n = B n , n = 1,...,N, that follows from ()2.33[) . implies the set of 
equations 

(0|$(0)|£ 2JV+ i- n ) = <0|$(0)|B n ) (4.22) 

for any local operator <3? of the massive minimal model. In the minimal models, and more 
generally in the sine-Gordon model, the matrix elements involving particles B n with n > 1 are 
related to the form factors (|3.2j) by residue equations on bound state poles of the type l)3.38j) . 
Due to ()2.33j) , in the minimal models these bound state equations constrain the form factors (|3,2j) 
themselves. For example, in the simplest case, N = 1, equation (|3.38f) holds with B 2 identified 
to B\ (see appendix E for the case of generic N). We find that the only way of satisfying (|4.22|) 
for any N when $ = TT is to take 13 

(0\K(0)\B n ) = (4.23) 

for n = 1, . . . , 2N , or equivalently 14 

lim r] n - x F^ [Q x + V ,e 2 + 2rj, . . . , 6 n + nrj) = (4.24) 

with Ok = 9\ + i(k — 1)£, k = 2, . . . , n and n > 1. 

The requirement of continuity in £ then leads to extend (|4.24|) to generic values of £. We 
checked explicitly up to n = 9 that the conditions (J4.4|) . (|4.6I) and (|4.24|) uniquely determine 
Fjf for generic ^. The explicit derivation up to n = 4 is given in appendix C, while appendix D 
contains the list of results up to n = 7. When inserted in ([4.12ft they provide the form factors 
of TT in the sinh-Gordon model seen as a perturbation of Liouville theory. The results apply 
to the sine-Gordon and the ^3-perturbed minimal models -M 2 /(2tv+3) as specified by i) and 
ii) above. Setting to zero the F^ T with n odd one obtains the result for the sinh-Gordon and 
sine-Gordon models seen as perturbations of the Gaussian fixed point, for which the reflection 
symmetry ip — > —ip holds. Few remarks about the free limit £ — > are contained in appendix F. 

5 Conclusion 

In this paper we identified the form factor solution corresponding to the operator TT in the 
sinh-Gordon model and in the ^i^-perturbed minimal models -M 2/ (2^+3)- The identification 
is obtained up to the arbitrary additive contribution of the operator dd@ which represents an 
intrinsic ambiguity in the definition of TT. 

We expect that the possibility of expressing the solution for TT as the superposition (|4.12f) 
of a three-particle kernel solution plus the contributions of the identity and the trace of the 
energy-momentum tensor together with its first two scalar derivatives is not specific to the class 
of models considered in this paper but is actually quite general. If so, the expressions of the 
coefficients a, d and e given here should be universal. 
13 For n = 1, 2 ijPSjl is implied by l|TTH for any value of £. 

14 The equivalence between 14.231 and (I4.24H for n > 2N follows from the periodicity at £ = £at. 
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The property, expressed by (|4.11l) and (|4.7j) . that TT behaves as a linear combination of 
derivative operators on states with more than four particles should also be very general. In 
particular, we remarked in [H] that it ensures that TT does not contribute to integrability 
breaking (at least to first order) when used to perturb a fixed point action. As a matter of 
fact, many examples are known of integrable massless flows in which TT is the leading operator 
driving the flow into the infrared fixed point |40| . 
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The local scalar operators of the theory satisfy H3.3j) - ()3.6j) with s$ = and S{6) = 1. The 
primary operators 



A Appendix 




(A.l) 




(A.2) 



n=0 



are subject also to the factorization property (|3.25j) . One then obtains 




(A.3) 



Ff = n\ (if) 



(A.4) 



Equation (|2.7j) gives in this free case 
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(A.5) 



and then 




(A.6) 




vrm 2 {Ff) 2 



(A.7) 



as the only non-zero form factors of this operator. Comparison with ()3.32j) fixes 




(A.8) 
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We write for ©. 
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The results (|'2.9|) and (|2.11j) then easily follow using (J4.9|) to evaluate over form factors the sum 
rules (|HIH1|) and [HI] 

X * = -i / # x (Q(%)HQ))conn ■ (A.9) 

B Appendix 

Here we derive the constraints imposed on the form factors of TT by the relations (|4.4|) and 
show how they lead to the property (|4.7j) . 

The r.h.s. of ()4.4|) is expanded introducing in between the two pairs of operators the res- 
olution of the identity (|4,9|) . Then the l.h.s and r.h.s. of 1)4.4(1 are rewritten in terms of 
form factors by iterative use of the crossing relation (|4.5|) . Let (m| = {B{6' m ) . . . B{9' 1 )\ and 
|n) = \B(8i) . . . B(6 n )) be the two states with the same energy and momentum in (|4.4j) . 

In order to avoid to sit directly on annihilation or bound state poles of the form factors we 
introduce an infinitesimal splitting, parametrized by r/, between the energies and momenta of 
these two states. The identity (|4.4(l is recovered in the limit 77 — > 0. 

The l.h.s. has an expansion of the form: 

(B(6' m + V ' m ) . . . B(9[ + 7701^(0)15(0! + m ) ■ ■ ■ B(6 n + Vn )) = 

m n 

F^ rn (0l ri +r]' m +i7r,..,e' 1 +r]' 1 +iTT,ei+r] 1 ,..,e n +r] n ) + 27r^^<5(?^ - rji + 6'j - 6i) 

j=i i=i 

3-1 i-1 

II II s(v'j -rf h + e- - e' h ) s(m - ^ + e k - e]) F^ m _ 2 (i,j) + ... (b.i) 

h=l k=l 

where rj' a = (n + a)i], a = 1, . . . ,m, 775 = 67/, b = 1, . . . , n, and F^ m _ 2 {h 3) is the form factor 
with n + m — 2 particles obtained omitting the particles B(9i + rji) and B{6'- + 77J). The dots 
in (|B.1|) represent terms which factorize p delta functions, a product of two-particle amplitudes, 
and form factors of TT with n + m — 2p particles with p = 2, Int((n + m) /2). 

Let us consider now the r.h.s. of (|4.4j) . Among the infinitely many terms generated by the 
insertion of Q4.9|) the only ones relevant for the identity (|4.4j) are those that do not depend on 
x for rj = 0. For each fixed k, {m\ T(x)\k) {k\ T(0)\n) — (m\ @(x)\k) (k\ B(0)|n) can be expanded 
in terms of the form factors by using expansions like ()B.1|) for each matrix element. Here the 
arguments of the delta functions are the differences between the rapidities of the states \m) or 
|n) and those of In these expansions the terms that do not depend on x are only those 
factorizing a set of delta functions which saturate all the integrations in (|4.9|) and fix the state 
\k) to one with the same energy and momentum of \n) and \m). Such terms are in finite number 
and are generated only if k = n or k = m, when the delta functions fix the state \k) to \n) or |m), 
respectively. We can now rearrange such terms according to their content in delta functions: we 
have a number of terms which do not factorize delta functions plus a number of terms which 
factorize one delta function, and so on as in formula ()B.1|) . The identity (|4.4j) is recovered 
imposing that in the l.h.s. and in the r.h.s. the terms that factorize the same delta functions 



18 



coincide in the limit rj — ► 0. An iterative structure appears for these identities (see (|4.16f) for the 
particular case m = n = 1). In general, the terms factorizing delta functions give the results of 
the identities 1|4.4[) with m — 1 and n — 1, while the terms without delta functions give the new 
condition 

min(n,m) 

limF^ m (9' m +rj' m +i7r,..,e[+7][+in,6i+ri 1 ,..,e n +T] n ) = lim{ ^ _ ^ ^ 

a 

II A (^ m - a+h - 0T n _ o+ J{[J ? £fn-2a(^ a +^_.+^^^ 
ft=l 



i ? ® +n _2a(^m-a+^-a+^'",^ 1 +^ 1 +^,6 l r 1 +??r 1 ,.-,6 l T„_ a +??r Il _ £i ; 



X 



^2®(C+C+^'-C- a+ i+i m -a +1 +^»^-a +1 +^„_ a+1 ,-^r n +r/rJ] + [T <-> T]} 



(B.2) 



where "P n is the group of permutations of n indices, r is a permutation and 

. . . | 1 for x = , 
A(x) = i 

[ for x ^ . 

Let us observe now that the r.h.s. of (|B,2|) goes always to zero at least as r/ 2 , because the p- 
particle form factors of T,T and O respectively factorize (cr^) 2 , (o£_i) 2 and x . The only 

exceptions to this situation arise if some rapidities are grouped in pairs that lie on annihilation 
poles or for rapidity configurations which intercept bound state poles. In every case, in terms 
of the parametrization (|3.7j) . the only way in which the form factors of TT can satisfy (|B.2|) is 
that for any p > 2 the function is the sum of terms each one factorizing f^i^l (&^2x) 
with some i and j such that i + j > 2. 

This property amounts to (|4.7jl provided that there are no cancellations with the denominator 
of 1)3. 7j) . This denominator can be written as 

/ 1 \ (p-i)/2 

[] COSh^-%) = [] (Xi+Xj) 

l<i<j<p \2 p cr p / l<i<i<p 

/ , \ (P"l)/2 



det L> (p) , (B.3) 



where Z)W is the (p — 1) X (p — 1) matrix with entries 



= ^ • (B-4) 



For p > 2, it is only for p = 4 that (jB.3|) is a sum of terms all containing non-vanishing powers 
of or and cancellations may occur. 
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C Appendix 



In this appendix we determine the form factors of the kernel /C up to four particles. At three 
particles the most general solution of the form factor equations satisfying 1)4.7)1 . 1)4.1 J) with I = 2 
and the initial condition 

F? = (C.l) 

16 



IS 



i5 + £$3 

Here, -F,^ 3 is defined by 



CJ3 CJ 3 



Qfa = (ffl<72 _ CT3 ) (C.3) 

in the parametrization 1)3.7)1 and generates the one dimensional space of the level solutions of 
the equations 1)3. 3)1 - 1)3. fij) which vanish at one particle. The first and second term in ()C.2j) are 
pure level 2 and level 1 parts, respectively, and 63, B 3 , C3 and D3 are dimensionless coefficients. 
The asymptotic factorization property fixes the level 2 solution, i.e. the coefficients 63, C3 and 
D 3 to 

63 = 6 = -(6) 2 , (C.4) 



C 3 = -6,D 3 = -6. (C.5) 
For n = 3 the condition 1)4.24)1 fixes the level 1 part, i.e. the coefficient B3 to 

B 3 = -6(1 + 2 cos 2^). (C.6) 

Now, let us consider the four-particle case. The most general level 2 solution to the form factor 
equations satisfying 1)4.7)1 and the initial conditions 

F$ = (C.7) 

is 

r 2^2 1 rp „2„ „ 1 p\ ^ ^2 1 n 



b 4 a 1 a 3 + E 4 a 1 a 2 (J 4 + D 4 <7 2 o" 3 + C 4 a 2 a 4 K aia 3 K K 

o + D 4 1± + A 4 1< A 

Here, F^ 4 is defined by 



4 --1--D4 -C4 --t-A 4 r 4 - . (C.8) 

a% o 4 



Q 4 4 = (crio-30-2 o-\a 4 ) 

and generates the one dimensional space of the level solutions of the equations 1)3.3)1 - 1)3.6)1 
which vanish at two particles. The asymptotic factorization property fixes the level 2 part of 
TT8|) . i.e. the coefficients b 4 to the same value of b, and C4, D4 and E4 to 

C 4 = -b, D 4 = 0, E 4 = . (C.9) 



3 We simplify the notation by dropping the superscript (n) on the symmetric polynomials. 
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The general result ()B.2|) for n = m = 2 is rewritten as 

lim Fj T {9 2 + r l + in,9 l + r l + in, 9 U 2 ) = lim{[F 2 T (^i + r] + in, 9 l )F^{9 2 + rj + in, 9 2 ) 

-F 2 S (9 1 + t] + in, 9 1 )F 2 9 (9 2 + r] + in, 9 2 ) - (Q)F®(9 2 + rj + in,9 x + V + in, e u 9 2 )] 

+ [T~T]}. (CIO) 

This expression implies for the four-particle kernel 

lim ff(0 a +r] + in,6 1 +ri + in,9 u 9 2 ) = 2[F 2 e (m, 0)] 2 (cosh2(^ - 9 2 ) - 1) , (C.ll) 

77 — >o 

where the r.h.s. comes from the identity 

F 2 T (9 1 +in,0 1 )F 2 T (9 2 +in,0 2 )-F 2 3 (9 1 +in,9 1 )F 2 9 (9 2 +in,9 2 ) = [ e 2 ^-^-l] [F®(in, 0)] 2 . (C.12) 

The above condition fixes the level part of (|C.8|) . i.e. the coefficient A4 to 

A 4 = 4&cos 2 £. (C.13) 

Finally, for n = 4 the condition ()4.24j) fixes the level 1 part, i.e. the coefficient B4 to 

B 4 = -26cos£(l + 4cos£ + 2cos2£). (C.14) 

This procedure has been implemented up to nine particles and the result is given in the next 
appendix. 

D Appendix 

We list in this appendix the functions which through (|H7j) . (|3.15j) and 

Q£ = -<e) 2 Q£ (D.l) 

determine in ()4. 12|) . The Qn have been determined explicitly up to n = 9. The functions 
Qf and Qg, however, are too cumbersome and we do not reproduce them here. We simplify 
the notation by dropping the superscript (n) on the symmetric polynomials. 

Qf = -2(Pi<4 - °l - ^3 - (1 + 2 cos[2£]) a 1 a 2 a z ){a x a 2 - a 3 ) (D.2) 
Qa = \{ol<?l - v 2 V4 ~ 2cos[£](l + 4cos[£] + 2cos[2^)cria 3 cJ4 + 4cos[£] 2 cr 2 ) 

CT 4 

(<7lV3V2 - of - £T 2 <T 4 ) (D 
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Qf = —{-ViVzvlvb + crjcr 5 (aja 3 - 4cos[£] 2 (7 3 £7 4 - a 2 a 5 ) + <7i 0-4(75 (- o\ 04 + (7 2 (2(2 + cos[2£])(7 2 

-(7 + 8cos[£] +8cos[2£] + 4cos[3£] + 2 cos[4£])ct 3 (75) + 4cos[£]((l + 4cos[£] + 2cos[2£] 

+ cos[3C])(7|(74 + (1 + 2cos[£] + 2 cos [2£] ) erf )) + erf (2(4 + 4cos[£] + 5cos[2£] + 2cos[3£] 

c 2 

+ cos[4^)a|(74(7| - (7 5 (crf(74 + 8cos[|] (2 + cos[£] + 2cos[2£] + cos[3£])(7 2 c7 5 + (7 3 (7 2 ) 
-(7 2 ((7l + 2(5 + 4cos[£] + 6cos[2£] + 2cos[3£] + cos [4£]) a z a\u h - afaf)) 
+(j?(2(7 4 (75((3 + 4cos[£] + 4cos[2£] + 2cos[3£] + cos[4^])<r| + (2 + cos[2£])(73(7 5 ) 

+ (7 2 ((73(73-4cOsK] 2 (7 5 3 ))) (D.4) 

Q6 = -o (^2(^3^5 ~ ^2^1 ~ 4C0S[^] 2 (7?(74(75 + 4cOS [£] 2 £7 2 (7f - £7 2 £7 2 C7 3 (7 5 £7 6 - £7 2 (7f(7 2 (7 6 
^6 

— 0"l(7 2 (74(7 2 (76 — £7^(73 £740"| (76 — 2cOs[£](7 + 16cOs[£] + 12cOs[2£] + 4cOs[3£] 

+ 2 COs[4^])(7 2 (7 2 (73£7 4 (7 2 (7 6 + £7^ £7 3 (7 4 (7 2 (7 6 + 2(6 + 11 COs[£] + 8 COs[2£] + 6 COs[3£] + 2 COs[4£] 

+ C0S[5C])(71(7|(74(7|C7 6 + 2(6 + 11C0S[£] + 8cOs[2£] + 6cOs[3£] + 2cOs[4£] + COs[5£])£7 3 (7f(7 2 £7 6 

+ 2(2 + COs[2^])(7i(7 2 (7|(7|£76 - 4cOs[£] 2 (7 3 (7 2 : £7 2 (76 - 4cOs[£] 2 (7f 2 Cf\ £7 6 + 2(6 + 11 COS [£] + 8cOs[2£] 

+6cos[3£] + 2cos[4£] + cos[5^])ala^a e + (5 + 4 cos[2^])(7?(73(7|(7 6 - 2(1 + 2cos[£]) 3 (l - cos[£] 
+ cos[2^])(7i(7 2 (73(7|(7 6 - 2(6 + 15cos[£] + 8cos[2£] + 8cos[3£] + 2cos[4£] + cos[5^])(7 2 (7 4 c7f a 6 
-(7 2 (7 4 (7f(7 6 + 4(2cos[£] + cos[3C])£7i(7|(7 6 + £7 2 <7 2 o- 3 (7 2 : (7 2 + 2(6 + 11 cos[£] + 8cos[2£] + 6cos[3£] 

+ 2cOs[4£] + C0S[5£])(7 2 (7 2 (7 2 £75(7 2 - 2(1 + 2 C Os[£]) 3 (l - COs[£] + C0S[2^])(71(7|(75(7 2 + 2(2 

+ COs[2£])(7 2 (7 2 a4(75(7 2 - 2(1 + 2cOs[£]) 3 (l - COs[£] + C0S[2^])(7 3 £7 3 (74(75(7| - 2<7i<720"30"4<7 5 <7g 

-(3 + 4cos[2£])(7 2 (7f(75(7 2 + 4cos[£] 2 a\ o\a\ — 2(6 + 15 cos[£] + 8 cos [2£] + 8 cos [3£] + 2 cos [4£] 
+ cos[5C])(7 3 (7 2 (7 2 (7 2 - (3 + 4cos[2^])c7icr|(7 2 (7 2 + 2(4 + 17cos[£] + 6cos[2£] + 9cos[3£] + 2cos[4£] 

+ C0S[5C])(7 2 (73(7 2 (7| + (1 + 2cOs[2£])(7 2 Cr 3 (7 2 (7 2 - 2(1 + CO S [2£] + COs[4£])<7i<74<7ij<7jj + (1 
+ 2 C0S[2^]) 2 (7|(7| - 4C0S[^] 2 (7 2 (7|(73(7| - £7 3 (7 2 (7 4 (7 3 ; + (1 + 2 C0s[2£])<jf <7 3 <7 4 <7jj + 4(2 COs[£] 

+ cos[3e])£7fc7 5 (7| - 2(1 + cos[2£] + cos[4£])(7 2 (7 2 (75(7 3 + 2(1 + 2cos[2£]) 2 (7l(7 3 (75(7g 

+(l + 2cos[2e]) 2 (7 3 £7^) (D.5) 
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Qj = — (afala 7 (4cos[^] 2 a 5 a 6 - (1 + 2 cos[2£])cr 4 cr 7 ) - crfa 6 (-a 7 (-a 3 a 4 a 5 aQ + (3 + 8cos[£] 

+4cos[2£] + 4cos[3£] + 2cos[4£])o- 3 , - 2(6 + 7 cos [2£] + 2 cos [4£]) a 5 a 6 a 7 

+8cos[£] 2 (l + 2cos[2^])cr 4 o-f) + a 2 a 6 (4cos[£] 2 a^a\ - cj 7 ((1 + 2 cos[2£])o- 4 cr 6 + 4cos[£] 2 o- 3 o- 7 ))) 

2 2^2 2 2 2 

-<7 6 c7 7 (4cos[£] a 2 a 3 a 7 -a 2 a 5 (a 3 a A + 4cos[£] cr 7 )-cr 5 ((l + 2 cos[2£]) a- 5 o- 6 -8cos[£] cos[2£]o- 4 o- 7 ) 

+4cos[C] 2 cr 3 (cr|cj5 + 4cos[£] 2 (l + 2cos[2^])o- 6 cj 7 ) + a 2 (-(l + 2 cos[2£])o- 3 cr 5 o- 6 + 4cos[£] 2 o- 2 

+<7 4 (of - 16cos[£] 4 cr 3 o- 7 ))) + cjicj 6 cr 7 (-cr|cr 4 cr 5 cj 6 - 4cos[£] 2 (-l + 2cos[£])(l + 2cos[^]) 3 ct|<j| 

+4cos[^] 2 cr|cr 6 cr 7 - 4 cos[£]of ct 4 (- (5 + 13cos[£] + 8cos[2£] + 6cos[3£] + 2cos[4£] + cos[5£])cr 5 o- 6 

+(1 + 2cos[£]) 3 (l - cos[£] + cos[2£])o- 4 cr 7 ) + (1 + 2 cos[2^])cr 3 (8cos[^ 2 cj| + (11 + 20cos[2£] 

+8cos[3£] + 10cos[4£] + 4cos[5£] + 2 cos[6£])cr 5 o- 6 cr 7 - 8cos[£] 2 (l + 3cos[2£] + 2cos[3£] 

+ cos[4£])o 4 o- 2 .) - cr|((3 + 4 cos [2£]) 0-5 ct 2 + cr 7 ((9 + 10cos[2£] + 2 cos[4£])cr 4 cr 6 - 8cos[^] 2 cr 3 cj 7 )) 

+ct 2 (2(2 + cos[2e])a|(J5(J 6 - 4cos[£] 2 (-l + 2cos[£])(l + 2 cos[^]) 3 cj 3 cj|cj 6 + (11 + 20cos[£] 

+20cos[2£] + 12cos[3£] + 10 cos [4£] + 4 cos [5f] + 2 cos [6£])o- 3 0^0-507 + 8cos[£] 2 (l + 4 cos[£] 

+3cos[2£] + 2cos[3£] + cos [4£]) of o- 6 o- 7 - <r 7 (-(19 + 26cos[2£] + 10cos[4£] + 2cos[6£])cr 2 

+4(3 + 5cos[£] + 4cos[2£] + 3cos[3£] + cos[4£] + cos [5£]) 050-7)) + 2(2(2 cos[£] + cos[3£])af <7 6 

-(3 + 3cos[2£] + cos[4^])o- 4 cr 5 cj| + 4cos[£] 2 o- 4 ofo- 7 + 2cos[£] 2 (l + 2 cos[2£])cr 7 (crf cr 6 + 2(1 

+2cos[£] + 2cos[2£] + 2cos[3^])<r?))) + of(2(13 + 18cos[«e] + 19cos[2£] + 10cos[3£] + 7cos[4£] 

+2cos[5£] + cos [6£]) erf 05 ct 2 ct 7 - 4cos[£] 2 (-l + 2cos[f])(l + 2 cos [£]) 3 o^oe of 

+o-|(-4cos[^] 2 cr 3 cr|cj 7 + (1 + 2 cos[2^]) 2 cr^) + <7 6 (-(l + 2cos[2£])of o- 6 of + cr 3 o- 5 cr 6 ((l 

+2cos[2£])of + 2(2 + cos[2£])o- 5 cr 7 ) + of(2(13 + 14cos[£] + 22cos[2£] + 12cos[3£] 

+13cos[4£] + 8cos[5£] + 5cos[6£] + 2cos[7£] + cos[8£])of + (19 + 26cos[2£] + 10cos[4£] 

+2cos[6^])a5CT 7 )) - 4cos[Z} 2 a 4 (alal + (1 + 2 cos[2£])<t 7 ((1 + 2 cos [£]) 2 o- 3 o| + 4cos[£] 2 of )) 

+o- 2 (4cos[£] 2 o| -2(7 + 22cos[£] + 11 cos [2^] + 12cos[3£] + 5cos[4£] + 2cos[5£] 

+ cos[6^])o- 5 cr|cj 7 - (3 + 4 cos[2£])ofo- 6 of + 2(7 + 22 cos[£] + 11 cos[2f] + 12 cos[3£] + 5 cos[4£] 

+2cos[5£] + cos[6£])o- 4 ofof + (1 + 2 cos[2^])cj 4 cr 5 cjf + o 3 a^{a A a^a\ - 2(3 + 3cos[2£] 

+ cos[4^])(jf))) + af{-a 3 al - 2cos[2f]<7 3 <7jj + 2ala 6 a 7 + 2 cos[2£]^ct 6 <t 7 - ofo- 5 o- 2 0- 7 

+Ua 3 a 5 ala 7 + 44 cos [£]o- 3 o- 5 o-fo- 7 + 22cos[2£]o- 3 <7 5 o|<7 7 + 24cos[3^]cr 3 cr 5 cr|cr 7 

+10 cos [4£]ct 3 0-506 a 7 + 4cos[5^]<j 3 a-5(j|o- 7 + 2 cos[6^]<T 3 <T 5 c7|cr 7 - 9a 3 a 2 a6cr 7 

-1Ocos[2£]<t 3 ct 2 0- 6 <7 2 - 2cos[4^]<j 3 <j|ct 6 (T 2 + 4(2cos[£] + cos [3£]) alow 3 - + 4of o- 6 o- 3 

+6cos[2^](7 2 c7 6 ^ + 2cos[4£]cj 2 cr 6 cr 7 i - 49a 2 a^ - 80cos[^]a 2 a^ - 80 cos [2£] of of 

-60 cos [3f] of of - 44 cos [4£] a 2 of - 32 cos[5£]ofof - 14 cos [6£] of of - 8cos[7£]ofof 

-2cos[8£]ofof - 20-50-7 - 2 cos [2^0-50-7 - ala 6 a 7 ((l + 2 cos[2£])of - 8cos[£] 2 (l + 4cos[£] 

+3cos[2£] + 2cos[3£] + cos[4£])ct 5 o- 7 ) + cr|(7 7 (2(13 + 18cos[^] + 19cos[2£] + 10cos[3^] + 7cos[4£] 

+2cos[5C] + cos[6^])o-|ct 2 - <7 5 ct 7 (4cos[£] 2 (-1 + 2cos[^])(l + 2 cos[£]) 3 ct 4 o- 6 + a 3 a 7 ) + 2a 2 ((2 

+ cos[2£])ct 4 o- 6 - (10 + 22cos[^] + 15cos[2^] + 12 cos [3f] + 6 cos [4^] +2cos[5C] + cos[6£])ct 3 ct 7 )) 

-4cos[^o- 4 ((6 + 10cos[C] + 10cos[2£] + 5cos[3^] + 2cos[4^] + cos[5C])o-|o-|a 7 + cos[^]a^ 
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+ (1 + 2cos[2e])(T 3 a|o-f - cos[e]o- 5 (o-^ + 4cos[e] 2 o- 3 o-?)) + <7 2 (2<7§<r 7 ((13 + 18cos[C] + 19cos[2e] 



+10cos[3£] + 7cos[4£] + 2cos[5£] + cos[6£])crf + (2 + cos[2£])<7 5 o- 6 cr 7 - 2cos[£] 2 ct 4 o- 7 ) + G-j{-a\o\ 
+(31 + 8cos[£] + 52 cos [2^] + 12 cos [3£] + 32 cos [4£] + 12 cos [5^] + 12 cos [6f] + 4cos[7£] 
+2cos[8^])<T4<t 6 <t? - 2((6 + 7 cos H + 2cos[4£])<x| + 8cos[£] 3 (2 cos[£] - 2cos[2£] 
+3cos[3£] + 2cos[4£] + cos[5£])ct 5 o- 2 ct 7 - 2cos[£] 2 cr 2 cr 2 )) + cr 3 (2(13 + 18cos[£] 
+19cos[2£] + 10cos[3£] + 7cos[4£] + 2cos[5£] + cos[6£])o-|o- 6 crf - o±(o\ + 2(18 + 22cos[£] 
+25cos[2£] + 12cos[3£] + 8cos[4£] + 2cos[5£] + cos[6£])ct 5 ct|ct 7 - trjjof) - <t 7 (ct|<7 6 
+8cos[£] 2 cos[2£] f r 7 ; ))))) (D.6) 



E Appendix 

In this appendix we explicitly derive the form factor equations which characterize the operator 
content of the $i ; 3-perturbed minimal model -^2,2^+3- The particles in such models obey the 
following bootstrap fusion algebra 

B a x B b — > B min ( a+b 2 N+i- a -b) and B a x B b ^ B\ a _ b \ , a, b = 1, . . . , N (E.l) 

which correspond to the bound state poles located at 



3 min(a+6,2iV+l-a-6) _ - (n , u \£n 

T 

2 



-i(a + b)^- (E.2) 
■ llh " -KT-\a-b\^), a^b. (E.3) 



The pattern (|E.1|) makes clear the possibility to describe all the particles and their fusions in 
terms of the fusion processes involving only the lightest particle B\. In this way all the bound 
state constraints can be handled inside the parametrization 1)3. 7 Jl . 

In particular, we can describe the particle ^niin(n,2Ar+i-n) by the fusion of n B\ particles set 
in pairs on the bound state pole of the B\B\ scattering channel. For the form factors of the 
generic local operator $ this implies 

(0|$(0)|5 min(n ,2iV+l-n)(£n)£i(0i) • • • B X (6' m )) = 



-i) 



n— 1 



lim lim e 1 --e n -iF m+n (9i + ei,..,9 n - 1 + 

€n— 1) 9 n , 9[,..,9' m ) (E.4) 

for 1 < n < 2N, 9 h = 9 X +i(h- i)£ N , h = 2, . . . ,n, and 6 n = 6% + i(n- l)&v/2. The result does 
not depend on the way in which the limit is done and T n can be written as 



T J UU^ 1 forl<n<iV, 
" ' Ut^U^-n^ 1 foriV<n<2iV, ' 



where T^ b is the three-particle coupling defined by 

Res fl=fl c S ab (9) =i(T c ab ) 2 . (E.6) 



24 



The use of formula (|E.4D allows to express the conditions Q4.22JI as 

lim... lim ei • ^AT-n^iV+l-n^l + e l) ••> #2iV-n + € 2N-n, O^N+X-n) = 

(-l) N ~ n „ Tn lim ... lim ei • ^l^Wi + ex, .., n _i + e n _ x , n ) , (E.7) 

J- 27V+l-n £ i^0 En-1— »0 

with ^ defined as above, 1 < n < N, 9' h = 9[ + i(h - 1)£ N , h = 2, . . . , 2N + 1 - n, and 
9[ = 9\ — i(2(N — n) + 1)^/2. Moreover, for each m + 2N B\ particles with m > 0, we can 
consider the case n = 2N in (|E.4|) which implies the equation 36 

lim ei • •e 2 Ar--iF* +27V (6' H — £x + ei,..,9 £n, 9i, .., 9 m ) = 

i( 2N -Vr 2N F* +1 (9,9 1 ,..,9 m ) . (E.8) 
In terms of the parametrization Q3.7J1 the above equations become 

Qm+l-ni^V ••' ^TV+l-n) = (h)Qn(6l, ■■, n ) , (E.9) 



Q%,+2n(0 H — 6v, • • • ! # + t;£n,Q - t;£,n, ■ ■ ■ ,0 — £n,0i, ■ ■ ■ ,0m) 



rN 

[I'll . ■ 

Qn {Ox,..., n ) 

and 

—$n,9 £m ..... 9 

2 2 2 

Vm (x,xx,... ,x m )Q m+1 (9,9x,... ,9 m ) , (E.ll) 

with 

N—X m N m 

V£(x, Xl ,...,x m ) = (-l) M ^ ±ii + 1 H [fcpaJW-l) "Q (x + Xi) "Q JJ (x _ _ XiW -fc } 

k=2 i=X k=2i=X 

(E.12) 

Xj = e ldi and w = e l ^ N . 

These bound state equations complete the equations (|3.18|) - (|3.21[) and characterize the oper- 
ator content of this model. Let us describe the case of the primary fields. For £n = 2tt/(2N + 1) 
a periodicity arises and the set of exponential operators e lk ^ v of the sine-Gordon model reduces 
to a finite set of 2N independent ones. Indeed, the form factors of e* fc ^ satisfy 

F^ 2N+1 H0i,..,0 n ) = F*{9 1 ,..,9 n ) , (E.13) 



F (2N+i)-k {e ^^ 9n) = ( _ 1} n+i F k^^ 9n) ^ (E . 14) 

implying that the independent among the operators e w can be chosen for k = l,..,2N. 
Now, only those with k = 1, ..,N have form factors satisfying the bound state equations (|E.9|I 
and (|E.11|) . Thus, the reduction of the operator content of the $i ; 3-perturbed minimal model 
M-2,2N+3 is implemented simply requiring the bound state equations. 

The form factors of the operator TT that we computed explicitly up to nine particles satisfy 
the bound state equations (|E.9|) and (|E. 1 1|) for any N when setting £ = as required by the 
property ii) in section 0J 
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F Appendix 



The expressions given in this paper for the form factors of the components of the energy- 
momentum tensor and for TT apply to the sinh-Gordon model seen as a perturbation of the 
Gaussian (C=l) fixed point when the form factors on states containing an odd number of 
particles are set to zero. In this case the limit 6 = corresponds to a free massive boson with 
background charge Q = 0. The form factors of and TT behave as 

Ff m ~ b 2 ^ , F% ~ b 2 ^ (F.l) 

as b — > 0. As expected, the only finite non-vanishing form factors at b = are obtained for 2 and 
4 particles, respectively; there are however infinities at particles for and at and 2 particles 
for TT that need to be subtracted. 

The usual finite trace operator with zero expectation value at b = is defined through the 
subtraction 1 

Q R = Q-(Q)I. (F.2) 
The divergences of TT are eliminated in the subtracted expression 

TT R = TT — am~ 2 d 2 B 2 Q — dm 2 — e m 4 1 + J 7 

= K + cddQ + r , (F.3) 

where a, d and e are the (infinite at b = 0) coefficients given in section 4, and T is a finite part 
needed to ensure that the asymptotic factorization (|4.6j) continue to hold for the regularized 

operators after the subtraction of the term proportional to d 2 B 2 Q. This implies that vanishes 
for n ^ 4 and is determined by 

2 2 

Of = <e> 2 ^ - 4 - a 2 a A ) , (F.4) 

so that 

Fl fR = {^f) 2 + 14^3 - 4a 4 ) . (F.5) 

If Fj TR is set to zero choosing c = 0, this is the only non- vanishing form factor of TTr. 

This result is easily compared with that for the matrix elements of T{x)T{x) — Q(x)Q(x) 
computed in free field theory with normal ordering regularization. The only difference arises 
in the coefficient of the term proportional to (Jia-s in l)F.5|) . amounting to a contribution of the 
operator ddfi 4 . This is precisely one of the derivative terms in (|4.3|) left unfixed by (|4.6|) , 



The limit b — » is understood in the equations below. 
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